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A direct spectroradiometric determina- 
tion of the temperature of freezing gold 
was performed by measuring the spec- 
tral radiances of a gold blackbody rela- 
tive to those of a laser-irradiated 
integrating sphere which was calibrated 
with absolute silicon detectors and an 
electrically calibrated radiometer. The 
measurements were performed at three 
laser wavelengths near 600 nm, and the 
temperature of the blackbody was calcu- 
lated by substituting the measured spec- 
tral radiances into Planck's radiation 
formula. The result obtained, 
7'au=(1337.33± 0.34) K, is 0.25 K be- 
low the gold-point assignment in the In- 
ternational Practical Temperature Scale 



of 1968 (IPTS-68) and has been adopted 
in September 1990 as the new gold- 
point value in the International Temper- 
ature Scale of 1990 (ITS-90). The effect 
of this change in the gold-point assign- 
ment on pyrometric, radiometric, and 
photometric measurement services pro- 
vided by the National Institute of Stan- 
dards and Technology is assessed. 
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1, Introduction 

In this paper, we describe a new spectroradio- 
metric measurement of the freezing temperature of 
gold based on absolute detector standards. Before 
describing our measurement, we review the defini- 
tions of this temperature in the International Practi- 
cal Temperature Scale of 1968 (IPTS-68, [1]) and 
in the International Temperature Scale of 1990 
(ITS-90, [2]), assess its significance for pyrometric, 
radiometric, and photometric measurement ser- 
vices provided by the National Institute of Stan- 
dards and Technology, and review gold-point 
measurements performed elsewhere since the adop- 
tion of the IPTS-68. 

The IPTS-68 allows the realization of tempera- 
tures, Tes, above the freezing temperature of gold 
by means of the equation: 



L^{T,,)/L^{T,,,^,)=:PQ.,T,,)/PQk,T,,,^;), 
where 



' 68,Au 



= 1337.58 K 



(1) 



(2) 



is the temperature assigned to the gold point; Lx{T) 
is the spectral radiance [Wcm"^sr~'] at wavelength 
A. of a blackbody of temperature T [K]; 



P(k,T)={c,/n ^?iMexp[c2/« XT]-\} 



(3) 



is the theoretical blackbody spectral radiance 
(Planck's radiation formula) for light with air 
wavelength X propagating in air with refractive 
index n ; and C\ and Cj are the first and second radi- 
ation constants. The gold-point temperature given 
in eq (2) was derived from gas thermometry. 
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The new ITS-90 was adopted by the Interna- 
tional Committee of Weights and Measures 
(CIPM) in September 1989 [2]. It supersedes the 
IPTS-68 and incorporates a new value for the 
freezing temperature of gold which is identical to 
the one reported in this paper [eq (35), below]. The 
ITS-90 extends the radiometric temperature defini- 
tion of eq (1) down to the freezing point of silver, 
and allows the realization of radiation temperatures 
in terms of any one of the silver [r9o.Ag= 1234.93 
K], the gold [r9o.Au= 1337.33 K], or the copper 
[79o.cu= 1357.77 K] freezing points. The T^o values 
of the silver, gold, and copper points are believed 
to be self-consistent to such a degree that the sub- 
stitution of any of them in place of the value of 
T'es.Au in eq (1) will not result in significant differ- 
ences in the measured values of T^q. 



2. Radiometric Significance of the Gold 
Point 

Realizations of the NIST scales of radiance tem- 
perature, spectral radiance, spectral irradiance, and 
photometric quantities, and the transfer of these 
scales to calibrated lamp standards, have been 
based on measurement chains derived from gold- 
point blackbodies for many years [3] ... [6]. Al- 
though the instrumentation and procedures used 
have undergone refinements in the meantime, the 
IPTS-68 gold-point assignment has remained the 
starting point for these scale realizations for the 
past 2 decades. Consequently, the change of the 
gold-point temperature will affect all of these 
scales. The magnitude of the scale changes may be 
estimated as follows. 

The realization of radiometric scales in terms of 
a primary gold-point standard involves measure- 
ments at specified wavelengths of the spectral radi- 
ance ratios of the gold-point blackbody and other 
blackbody sources which may then be used as sec- 
ondary standards over wide temperature and 
wavelength ranges. The basic measurement equa- 
tion for these radiance transfers is eq (1), which 
allows the determination of the temperature of a 
secondary standard source from the measured radi- 
ance ratio LxiT)/LJT/^^ and the assigned value of 
Tau. In the Wien approximation of Planck's radia- 
tion formula. 



P(k,T)~(cx/n^\^) exp(-C2/«Xr), 



we obtain 



(4) 



L,(r)/L,(rAu)~ exp[feA)(l/rA.- \/T)] (5) 

and, hence, by implicit differentiation with respect 
to T and Tam, 



dT/T'-dTj,yTi,~0, 



(6) 



if the radiance ratio on the left-hand side of eq (5) 
has a fixed value. This shows that a radiance tem- 
perature T based on a gold-point value Tau will be 
replaced by T + AT, where 



Ar~(r/rAu)'ArAu, 



(7) 



on a redefined scale in which T^^ is replaced by 

When this result is substituted into the Wien ap- 
proximation for the spectral radiance of the black- 
body at temperature T, we obtain 

AP(X,T)=[dP(X,T)/dT]AT~P(X,T)c2AT/(XT^) 

~PQ,,T)c2AnjQ,TU. (8) 

Thus, the relative change in spectral radiance due 
to a change AT^au in the gold point is 

AP(\,Typ(\,T)~c2 ArAuAxrij 

= IA4X10-' ATUQ^TL), (9) 

where T is expressed in kelvins and X in meters. 
This result is independent of the temperature of the 
blackbody and holds for spectral irradiance mea- 
surements as well, since radiance/irradiance 
transfers are based on purely geometrical consider- 
ations. 

The NIST photometric scales (luminous inten- 
sity and luminous flux) are realized and maintained 
at r=2856 K (CIE Source A) by direct compari- 
son with the NIST spectral irradiance scale. The 
effect of a reassignment of the gold point may be 
estimated by evaluating eq (9) at the Crova wave- 
length for Source A (X=570 nm) [7]. This gives 



Aev/ev= 1.41 X 10-2 AT^Au, 



(10) 



where ATau is again expressed in kelvins, and Q^ 
denotes a generalized photometric quantity. 

From eqs (7), (9), and (10), it can now be esti- 
mated that the -0.25 K change in the IPTS-68/ 
ITS-90 gold-point assignments will result in the 
approximate changes of reported calibration values 
listed in table 1. Although these changes are well 
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Table 1. Quoted uncertainties (3or) of current NIST measure- 
ment services for radiation thermometry, radiometry, and pho- 
tometry and estimated changes in reported values resulting from 
a —0.25 K reassignment of the gold-point temperature 



Quantity 


Uncertainty 


Change of 
value 


Radiance temperature [4] 






800 °C 


±0.5 °C 


-0.2 °C 


1100 


0.6 


-0.2 


1400 


0.8 


-0.4 


1800 


1.3 


-0.6 


2300 


2.0 


-0.9 


Spectral radiance [3] 






250 nm 


±1.6% 


-0.8% 


350 


1.2 


-0.6 


650 


0.6 


-0.3 


900 


0.5 


-0.2 


1600 


0.4 


-0.1 


2400 


0.4 


-0.1 


Spectral irradiance [5] 






250 nm 


±2.3% 


-0.8% 


350 


1.4 


-0.6 


650 


1.0 


-0.3 


900 


1.3 


-0.2 


1600 


1.9 


-0.1 


2400 


6.5 


-0.1 


Luminous intensity [6] 






2856 K 


±1.0% 


-0.4% 


Luminous flux [6] 






2856 K 


±1.4% 


-0.4% 



within the quoted uncertainties for routine NIST 
measurement services, they will have an effect to 
reconcile small discrepancies between radiometric 
scales that have been observed in recent years. For 
example: In an intercomparison between indepen- 
dent radiometric scales based on silicon-photodi- 
ode physics, gold-point blackbody radiation, and 
synchrotron radiation, Schaefer, Saunders, and 
Hughey [8] found that the blackbody scale indi- 
cated a spectral irradiance at 600 nm which was 
about 0.8% higher than the detector scale, and 
about 1% higher than the synchrotron scale. The 
— 0.25 K adjustment of the gold-point temperature 
improves the overall agreement of these intercom- 
parisons. 

In an international intercomparison of photomet- 
ric base units, conducted by the Comite Consultatif 
de Photometrie et Radiometrie (CCPR) in 1985 
[9,10], the NIST luminous intensity data were the 
only ones derived from a gold-point based mea- 
surement chain. They fell within the spread of the 



intercomparison, but were 0.5% higher than the 
average of the data reported by 14 other national 
laboratories, all of which had realized the candela 
with electrically calibrated radiometers. It may be 
seen from table 1 that the NIST data would have 
been within 0.1% of the world mean, had the gold- 
point temperature been 0.25 K lower. 



3. Previous Radiometric Determinations 
of the Gold Point 

The first measurement which projected a differ- 
ence between the thermodynamic temperature of 
the gold point and its assigned IPTS-68 value was a 
precision determination of the Stefan-Boltzmann 
constant which was published by Blevin and 
Brown [1 1] of the Australian National Standards 
Laboratory (NSL) in 1971. Based on the IPTS-68 
gold-point assignment their experiment, in which 
an electrically calibrated radiometer (ECR) was 
used to measure the spectrally total radiance of a 
blackbody radiator at the freezing temperature of 
gold, yielded a result which was about 0.1% lower 
than the value of the Stefan-Boltzmann constant 
(Iv^/ISh^c^) calculated from the best values of the 
speed of light, c, and the Planck and Boltzmann 
constants, h and k. Blevin and Brown suggested 
that their experiment could be regarded as a radio- 
metric determination of the gold point based on the 
theoretical value of the Stefan-Boltzmann constant, 
the corresponding result for the gold point being 
1337.27 K, about 0.3 K lower than the IPTS-68 
value. 

Guildner and Edsinger [12] of the U.S. National 
Bureau of Standards noted that this result was con- 
sistent with their own work, which had revealed 
that sorption errors in gas thermometry caused a 
systematic difference between thermodynamic and 
IPTS-68 temperatures. Assuming that the effects of 
sorption were approximately linear with wave- 
length, they predicted that the temperature of 
freezing gold on the thermodynamic scale will 
prove to be lower than the value on the IPTS-68 
by at least some 0.1 K. They suggested, however, 
that the IPTS-68 should not be changed without 
better understanding of the thermodynamic scale, 
as realized experimentally. 

In 1975, Bonhoure [13] of the Bureau Interna- 
tional des Poids et Mesures (BIPM) in France 
reported pyrometric measurements of the tempera- 
ture differences of variable-temperature blackbod- 
ies in the 904 to 1338 K range, one kept at the 
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IPTS-68 antimony temperature (903.89 K) and the 
other at selected higher temperatures. Assuming 
the antimony point to be correct, Bonhoure found 
that his measurements were in 0.05 K agreement 
with the IPTS-68 near the gold-point. Andrews 
and Gu [14] repeated this measurement at the Na- 
tional Physical Laboratory (NPL) of the U.K., us- 
ing a photon-counting optical pyrometer to 
compare a gold-point blackbody to a reference 
source at the antimony point. They, too, concluded 
that the IPTS-68 gold-point value was not in error. 

In the meantime, Guildner and Edsinger [15] had 
achieved significant improvements of high-temper- 
ature gas thermometry and showed that the ther- 
modynamic temperatures from the triple point of 
water up to 730 K were consistently lower than 
the IPTS-68 values. Jung [16] of the German 
Physikalisch-Technische Bundesanstalt (PTB) car- 
ried this work further by performing pyrometric 
measurements at the antimony point relative to a 
blackbody at a thermodynamic temperature of 
729.15 K, as established by gas thermometry. His 
work showed that the IPTS-68 antimony-point as- 
signment was in fact too high by 0.15 K, thus ne- 
cessitating numerical corrections to the gold-point 
values obtained by Bonhoure and by Andrews and 
Gu. The corrected values, as reported by Tischler 
and Rebagliati [17] and Hudson et al. [2], respec- 
tively, were 1337.20 and 1337.33 K. An expression 
equivalent to eq (7), above, was used to make the 
corrections. 

Tischler and Rebagliati [17] also recalculated the 
results of five other researchers who had per- 
formed pyrometric measurements of the thermody- 
namic temperature difference between the silver 
and gold points between 1973 and 1982. To allow a 
uniform analysis, Tischler and Rebagliati inter- 
preted all of these measurements as being per- 
formed with the silver point as the reference. The 
five gold-point values so recalculated ranged from 
1337.18 to 1337.29 K, their average being 1337.22 
K. The silver-point reference temperature assumed 
was 1234.88 K, a value which is only 0.01 K lower 
than more recent results obtained by Jones and 
Tapping [18] and Fischer and Jung [19], and 0.05 K 
lower then the ITS-90 value[2]. 

The data described above are summarized in 
table 2, with uncertainties quoted at the three stan- 
dard-deviation level. The values given are all 
lower than the IPTS-68 gold-point assignment. 

It should be noted, however, that the only direct 
determination of the gold point is that by Blevin 
and Brown [11], which involved absolute radio- 
metric measurements and interpretation of the data 



Table 2. Results and uncertainties (3cr) of gold-point measure- 
ments performed since 1971 



Authors 



Reference 
Temperature 



Result (K) 



Blevin and Brown [11] None 1337.27±0.4 

Bonhoure [13], [17] Antimony 1337.20±0.5 

Andrews and Gu [14], [2] Antimony 1337.33±0.3 

Tischler and Rebagliati [17] Antimony/Silver 1337.22±0.5 



obtained in terms of a fundamental law of physics. 
The others are pyrometric measurements of radi- 
ance ratios which relied, directly or indirectly, on 
the IPTS-68 assignment of the antimony point and 
agreed with the Blevin and Brown measurement 
only after allowance was made for an error in the 
antimony assignment. We conclude that radiomet- 
ric and pyrometric temperature measurements at 
higher temperatures are, in themselves, reliable but 
may be encumbered with uncertainties if they are 
based on lengthy measurement chains which in- 
volve lower-temperature reference points. Our 
own measurement is similar to the Blevin and 
Brown experiment in that it is a direct measure- 
ment based on absolute radiometry, but differs 
from it in that we measured the spectral radiance of 
a gold-point blackbody and used Planck's law to 
derive its temperature. 



4. Apparatus 

4.1 Overview 

Our technique consisted of measuring the spec- 
tral radiance of a uniform, monochromatic source 
using two independent radiometric scales estab- 
lished by a standard silicon photodiode and an 
ECR. As the uniform, monochromatic source we 
used the exit port of an integrating sphere irradi- 
ated by monochromatic laser light near 600 nm. 
Then, by comparing the spectral radiance of this 
source to that of a gold-point blackbody with a 
monochromator-based spectroradiometer, we de- 
termined the blackbody radiance and calculated 
the gold-point temperature from Planck's radiation 
formula. The use of independent detector standards 
provided redundancy of the measured data and in- 
creased the likelihood of high accuracy and inter- 
nal consistency of the results. 

Figure 1 shows a diagram of the arrangement of 
the equipment used in these measurements. The 
three detector systems (silicon diode, ECR, and 
spectroradiometer) are mounted on a movable 
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Gold 



Insulation 



Wavemeter 
Dye Laser 



Ar, Kr, or HeNe 
Laser 



Figure 1. Gold-point measurement apparatus. 

carriage which permits each of them to be moved 
into position to view the laser-irradiated integrat- 
ing sphere. In addition, the carriage can be posi- 
tioned so that the spectroradiometer views the 
gold-point blackbody. 

Also indicated is a silicon monitor diode which is 
fixed above the sources and detectors and views 
the laser sphere from slightly outside the field of 
view of any of the other detectors. The purpose of 
this monitor was to compensate for small fluctua- 
tions of the laser-sphere source and for the pres- 
ence of a neutral-density filter which was inserted 
in the beam for the silicon-diode and spectrora- 
diometer measurements, as described in section 4.7. 

Most of the simple operations of the equipment 
are computer controlled. This includes controlling 
the blackbody heater current, positioning the mov- 
able carriage, operating the spectroradiometer 
wavelength drive, triggering the simultaneous mea- 
surements by the monitor detector and the other 
detector systems, operating the shutter for dark- 
current measurements, and setting the integration 
time for the DVM readings (usually 5 s). 

Throughout the discussion which follows, we 
state numerical estimates of potential sources of 
systematic error. These estimates are upper limits 
of systematic errors and will be combined in 
quadrature with the three-sigma random uncertain- 
ties of our measurements. 

4.2 Gold-Point Blackbody 

The blackbody construction is sketched in fig- 
ures 2 and 3. The furnace incorporates a cylindrical 
heat pipe between the electrical heater and the alu- 
mina tube which holds the crucible of gold. The 




Water Cooling 
Figure 2. Cross section of heat-pipe blackbody furnace. 




8 mm 



3 mm. dia. 



Figure 3. Blackbody inner cavity dimensions. 

crucible and inner cavity are made of graphite, as 
are many of the rings which limit the viewing cone 
in front of the crucible and support the thermocou- 
ple tube behind the crucible. A slow flow of argon 
is directed from the back to the front of the furnace 
within the alumina tube, and fittings on its ends 
permit maintaining a slight overpressure of argon 
except when the front shutter is opened for radi- 
ance measurements. The front graphite rings are 
regarded as sacrificial, and after five or ten melt- 
freeze cycles a few of the outermost rings had gen- 
erally become partially consumed by oxidation and 
were replaced. 

A microcomputer controls the heater power sup- 
ply by turning the furnace on, ramping the power 
up gradually over several hours, and finally alter- 
nating power levels to produce melts and freezes 
which typically last about 45 min each. The ther- 
mocouple provides the computer with the neces- 
sary continuous information about the state of the 
furnace which governs the succession of power 
changes. Spectroradiometric data from a represen- 
tative melt-freeze cycle are shown in figure 4. 

To confirm the radiometric adequacy of the 
blackbody design, several of the easily changed 
design parameters were tested by performing spec- 
troradiometric matches between an ultra- stable 
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D 50 100 150 

Time [minutes] 

Figure 4. Spectroradiometric blackbody melt and freeze curve 
(514 nm). 

vacuum tungsten strip lamp and the blackbody 
melts and freezes with modified blackbodies: 

(1) Furnaces with different internal diameters were 
compared. The two designs provided// 10 and// 
6.5 viewing solid angles and crucibles containing 
0.7 and 1.3 kg of gold, respectively. 

(2) Two different samples of nominally 0.99999 
pure gold, acquired at different times, were com- 
pared in the smaller furnace, and a third sample in 
the larger furnace. 

(3) Three different inner blackbody cavity open- 
ings were compared. Their diameters were 1.0, 1.5, 
and 2.0 mm. 

None of these changes affected the spectral radi- 
ance of the cavity within 0.01%. For the measure- 
ments reported here, the larger furnace with the 
crucible containing 1.3 kg of gold and 1.5-mm di- 
ameter aperture was used. For the calculation of 
the gold-point temperature from our spectral radi- 
ance measurements and the Planck function, we 
used an emissivity of 0.9999+0.0001, calculated 
from the dimensions of the cavity by Quinn's and 
Ford's technique [20,21], assuming uniform tem- 
perature. 

Since heat flows from the outer heat pipe to the 
inner cavity during melts, and tends to flow from 
the inside to the outside during freezes, we ex- 
pected differences in the relative spatial distribu- 
tions of the blackbody radiance during melts and 
freezes, as well as changing distributions during the 
course of a melt or freeze. Relative radiance distri- 
butions were measured during melts and freezes 
carried out for this purpose, by moving the spec- 
troradiometer in 0.6-mm steps across the face of 



the blackbody furnace, and by raising and lowering 
the furnace in 1.8-mm steps. Generally, four com- 
plete mappings could be acquired during the 
course of a single melt or freeze. We show in figure 
5 examples corresponding to different times during 
a melt or freeze. In order to eliminate the effect of 
the change in the distribution during the 10 or 15 
min required for each mapping, at each point of the 
source plane the four spectroradiometer readings 
obtained during the complete melt or freeze were 
fitted to a linear function of the time of measure- 
ment. Then, "instantaneous" mappings were con- 
structed by interpolation and used in the 
computation of the spectroradiometer spatial ap- 
paratus function correction /h defined in section 
5.2.3. 

Blackbody Aperture 




-12 -10 -8-6-4-2 2 4 6 8 10 12 
Scan Position [mm] 

Figure 5. Relative spatial distribution of blackbody radiance 
(647 nm) at different times during melts and freezes, (a) end of 
melt (b) beginning and middle of melt (c) beginning of freeze (d) 
end of freeze. 



4.3 Laser Sphere 

The laser-sphere source is a 50-mm diameter 
copper integrating sphere with a 6-mm diameter 
entrance port located about 120° from a 20-mm di- 
ameter exit port. The interior of the sphere is 
coated with a 3-mm thick layer of compacted 
Halon (polytetrafluoroethylene) powder [22]. On 
the outside of the sphere, a single loop of copper 
cooling water tubing is soldered to each hemi- 
sphere. 

A 20-W argon ion laser, a 5-mW helium-neon 
laser, and a 6-W krypton ion laser were used to 
illuminate the sphere at 514, 633, and 647 nm. The 
air wavelengths assigned to these laser transitions 
in Ref. [23] are given in table 3. The power output 
of the lasers remained stable to better than 1% over 
several hours, with no further stabilization than the 
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Table 3. Laser wavelength used in this experiment [19] 



Laser 



Transition 



Air wavelength 
(fj,m) 



Ar II 3/)*(Jp)4;pC'D?/2)->3;p''('P)4ifP3/2) 0.514533 

HeNe 2;p'CPi/2)5i[l/2]?-^2;p5(2Pi/2)3/)[3/2]i 0.63281646 
Kr II 4/)''('P)5/)(*P5/2)->4/)''('P)5i(3p3/2) 0.647100 



built-in standard power monitoring system of the 
laser. Since all of our measurements on the laser- 
sphere were made relative to a simultaneously 
recorded signal from the monitor detector, this sta- 
bility was adequate. The helium-neon laser pro- 
vided too little power for accurate ECR readings, 
so at 633 nm we only compared silicon-detector 
and spectroradiometer measurements. The laser 
beams were directed by a series of plane mirrors 
through a quarter-wave plate to a lens which fo- 
cused the beam to a point in the sphere entrance 
port. 

Because the spectroradiometer is sensitive to lin- 
ear polarization, it was necessary to know the po- 
larization of all sources. The blackbody can safely 
be assumed to be unpolarized, but we felt less cer- 
tain about the laser-sphere source since the laser 
beam is so highly linearly polarized. Consequently, 
we inserted the quarter-wave plate in the laser 
beam to convert its linear polarization to circular 
polarization entering the sphere. Then we mea- 
sured the linear polarization of the laser-sphere 
source by the use of a rotating sheet polarizer and, 
as a polarization indifferent detector, another inte- 
grating sphere with a silicon detector in its exit 
port. With careful attention to alignment, this tech- 
nique yielded consistent values of 0.0007-1-0.0003 
for the degree of linear polarization of the sphere 
output. This is sufficiently small, so that we treated 
the laser-sphere source as unpolarized. We in- 
cluded in our error analysis a polarization uncer- 
tainty of ±0.02%, which is the product of the 
measured degree of polarization of the sphere out- 
put and the polarizance of the spectroradiometer 
(0.3) in the spectral region of interest. 

The relative spatial radiance distribution of the 
laser-sphere source was measured with the spectro- 
radiometer. The exit aperture of the sphere was im- 
aged with a tiltable spherical mirror into the source 
focal plane of the spectroradiometer. Then, the 
horizontal scanning capability of the spectrora- 
diometer, coupled with the vertical tilting motion 
provided by the mirror, permitted a complete map- 



ping of the 255-mm^ source area in a 1 X 1 mm grid. 
Initial attempts to perform this mapping showed a 
fine structure which depended strongly upon laser 
beam positioning in the entrance port. This proved 
to be caused by a fine speckle pattern in the sphere 
output which we were able to smooth out by rotat- 
ing one of the plane mirrors to introduce a small, 
rapid wobble of the beam during the measure- 
ments. The spatial distributions thus obtained are 
reproducible, smooth, and independent of laser 
alignment and wavelength. In figure 6, we show 
the distribution at 514 nm. 

Since the spectroradiometer views the laser 
sphere within a larger solid angle (approximately 
//12) than either the silicon diode {f/\9) or the 
ECR (f/25), we required that the laser-sphere 
source be Lambertian. We tested for uniformity of 
the irradiance field in the detector plane by trans- 
lating the silicon-diode detector horizontally and 
vertically in the detector plane and recording the 
ratio of the signal obtained to that of the monitor 
detector. The cos'*0 dependence expected for a 
Lambertian source was obtained to within -1-0.01% 
over our largest (f/^2) viewing angle. 

4.4 Spectroradiometer 

The spectroradiometer used in our measure- 
ments was a Cary-14 prism-grating double 
monochromator'. The fore optics indicated in 
figure 1 consisted of a plane mirror and a 1-m 
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Figure 6. Relative horizontal distribution of laser sphere radi- 
ance at exit aperture (514 nm). 



' Certain commercial equipment, instruments, or materials are 
identified in this paper to specify adequately the experimental 
procedure. Such identification does not imply recommendation 
or endorsement by the National Institute of Standards and 
Technology, nor does it imply that the materials or equipment 
identified are necessarily the best available for the purpose. 
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radius-of-curvature spherical mirror. These were 
positioned to produce 1:1 images of sources on the 
entrance siit, while keeping the optical axis within 
less than 3° of the axis of the spherical mirror. A 
circular aperture stop at the spherical mirror lim- 
ited the angular field of view of the radiometer to 
f/Yl (except in the horizontal direction, in which it 
was limited to //1 6 by the internal optics of the 
monochromator). A polished metal mask contain- 
ing a 1.0-mm wide by 0.8-mm high rectangular 
opening was moiinted just in front of the entrance 
slit. Since the slit width used for all measurements 
was 0.6 mm, the effective field stop of the spectro- 
radiometer system was a rectangle 0.6-mm wide by 
0.8-mm high. A telescope focused on the entrance 
slit mask gave the operator a magnified view of the 
source image at the slit, to help in the positioning 
and focusing of sources. In the detector compart- 
ment behind the exit slit of the monochromator, a 
photomultiplier and a small HeNe laser were 
mounted on an externally positionable sliding plat- 
form, so that either can be placed on axis. The 
photomultipher is a thermoelectrically cooled, 
magnetically shielded, end-on, eleven-stage EMI 
9658 tube with S-20 response. The laser is used to 
define the optic axis of the spectroradiometer sys- 
tem for the positioning and orienting of sources. 

Measurements of the spectroradiometer slit-scat- 
tering function have been described previously 
[24], They employed a laser-sphere source and con- 
sisted of measurements over a broad range of dye 
and ion laser wavelengths while the monochroma- 
tor remained set at a fixed wavelength. These mea- 
surements were supplemented by information 
obtained by scanning the monochromator wave- 
length setting while viewing a constant ion laser 
line. An example of the resultant slit-scattering 
function is shown in figure 7. Using these slit-scat- 
tering functions and the definition of eq (31), be- 
low, we can calculate the spectral slit width. 
During the course of the measurements reported 
here, the slit-scattering function was measured at 
each wavelength by performing abbreviated 
monochromator scans which covered only about 
10 nm. Consequently, the spectral slit widths calcu- 
lated from these data were corrected for the trun- 
cation of the far wings by multiplying by a factor 
1.0025+0.0004, obtained from the earlier, com- 
plete scans by similarly truncating those data. For 
our nominally 0.6-mm slit widths, we obtained 
bandwidths in the neighborhood of 2.5 nm. The 
factor cr(X,XL) in eq (31) causes the slit width to 
depend on the precise coincidence of the 



monochromator wavelength, X, at the laser wave- 
length, Xl- The uncertainty of this coincidence 
(+0.05 nm) is included in the composite random 
uncertainty of our measurements, as discussed in 
section 6.6. 

To measure the spatial apparatus function of the 
spectroradiometer, we aimed it at the center of the 
laser-sphere source and mounted a mask containing 
a 0.3-mm diameter pinhole in the focal plane just in 
front of that source. Then we recorded the spectro- 
radiometer signal as the mask was translated hori- 
zontally and vertically across the source aperture. 
These traverses were then separately normalized to 
unity at the center. Examples are shown in figure 8. 
We found the spatial apparatus function to be es- 
sentially independent of wavelength over the spec- 
tral range used in our experiment. 
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Figure 7. Slit-scattering function of spectroradiometer (647 
nm). 
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Figure 8. Spatial apparatus function of spectroradiometer (514 
nm). 
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Although our gold-point measurement does not 
rely on an absolute calibration of the spectrora- 
diometer, the spectral responsivity Rs(X) of the in- 
strument is needed to evaluate the slit-function 
correction of the measured blackbody radiance [eq 
(32), below]. We determined RsQ^) with an accu- 
racy sufficient for this purpose by calibrating the 
spectroradiometer with a standard tungsten-strip 
lamp. The result of this calibration is shown in fig- 
ure 9. 
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Figure 9. Relative spectral response of spectroradiometer. 

4.5 Silicon Diode 

The silicon diode detector is a windowless, 1- 
cm^ p-n type, Hamamatsu 1337-lOlOQ diode oper- 
ated without temperature control. The response of 
this detector type has been found to be quite inde- 
pendent of temperature at ambient temperatures in 
the 500 to 700 nm spectral region [25]. The spatial 
uniformity of the diode response was mapped over 
a 7x7 mm grid with a 633-nm laser beam and was 
found to be better than ±0.1%, as shown in figure 
10. 

The silicon diode was calibrated by senior staff 
of the NIST Photodetector Metrology Project 
with a stated three-standard-deviation uncertainty 
of ±0.2%. This calibration was performed using a 
monochromator-based spectral response compara- 
tor and a high-accuracy, laser-based detector 
characterization facility. The diode spectral re- 
sponsivities at 514, 633, and 647 nm were measured 
by comparison to three fully characterized, abso- 
lute QED-200 detector standards. The QED stan- 
dard is a light-trap configuration of n-p type, UDT 
UV-100 silicon diodes which, by significantly re- 
ducing reflection losses, permits a simplified, accu- 
rate self-calibration procedure [26]. During our 
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Figure 10. Relative horizontal distribution of silicon-diode re- 
sponsivity (633 nm). 

measurements, we observed no changes in the 
diode response with respect to the ECR and moni- 
tor-diode readings. The overall agreement of the 
silicon-diode and ECR readings taken at 514 and 
647 nm was on the order of 0.05% (see table 7). 



900 4.6 Electrically-Calibrated Radiometer 



The ECR used was built by the National Insti- 
tute of Measurement and Testing Technology 
(NIMTT) of the P.R.C. and is similar to that de- 
scribed in reference [27]. As shown in figures 11 
and 12, it has two identical conical copper receiv- 
ers of 11° 19' half angle, 25-mm length, and 10-mm 
opening. The heater windings are made of 0.03 mm 
diameter manganin wire, are uniformly and tightly 
glued to the insides of the cone, and are coated 
with black paint and kerosene smoke. The outside 
of each cone is equipped with a uniform, circular 
array of 72 nickel-chromium, constantan (type 
"E") thermocouples made of 0.1 -mm wire. Only 
one cone is exposed to radiant or electrical heating, 
respectively, and the thermocouple arrays on the 
two cones are connected oppositely. In this man- 
ner, only the temperature difference of the cones is 
registered and environmental temperature changes 
(such as those caused by adiabatic air pressure fluc- 
tuations [28]) are largely canceled. The laser- 
sphere aperture and its distance from the ECR 
were chosen so that the irradiated and electrically 
heated portions of the ECR cone were closely 
matched. 

Due to residual differences in the electrical and 
radiant heating of an ECR, the temperature rise 
produced by radiant power is not exactly the same 
as that produced by an equal amount of electrical 
power. This non-equivalence of electrical and 



57 



Volume 95, Number 1, January-February 1990 

Journal of Research of the National Institute of Standards and Technology 



Heater 




Thermocouples 
in series 



Figure 11. Cross section of ECR. 
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Figure 12. ECR cone details. 

radiant heating has been examined in detail by 
Hengstberger [29]. The determination of the small 
correction for this non-equivalence is generally a 
one-time characterization of the ECR, as con- 
trasted with the routine calibration by electrical 
substitution heating to achieve a thermocouple sig- 
nal which matches that of an unknown radiant in- 
put power. Our characterization of the ECR was 
based on the equations 



'ECR- 



-Recr P'> 



P'=P(l-|-K-p-T7-e), 



(11) 
(12) 



where Secr is the ECR signal (thermocouple 
EMF), i?ECR is its responsivity, P' is the thermal 
power at the thermocouple hot junctions, P is the 
total (radiant or electrical) power applied to the 
ECR, and the remaining quantities are small cor- 
rection terms: 

K is the environmental heating correction, 
p is the reflection/emission loss, 
Tj is the air convection/conduction loss, 
€ is the sum of other losses. 

The last term, e, denotes the sum of all losses (such 
as thermal emission from the back side of the cone) 
which are identical for radiant and electrical heat- 
ing and, therefore, do not require a correction. The 
other correction terms depend on the type of heat- 
ing and will accordingly be distinguished by sub- 
scripts, r and e, so that k^ denotes radiant case 
heating, k^ denotes electrical lead heating, and so 
on. The characterization of the instrument consists 
of measuring the six correction terms k„ k^, p„ Pe, 
T/r, and Tje one at a time in specially designed exper- 
iments. To measure k„ we covered the ECR aper- 
ture with a styrofoam plug and measured the 
thermopile signals with and without a radiant 
power Pr falling onto the face of the instrument. 
The lead heating effect was determined by rear- 
ranging the ECR circuitry so that the heat gener- 
ated by the leads alone could be measured in 
addition to the heat generated by the heater wind- 
ings and leads. The reflection/emission loss was 
measured by irradiating the ECR with a laser beam 
entering through the small central hole in a spheri- 
cal mirror which could be rotated to focus the radi- 
ation remitted by the ECR back into the receiver 
cone, or to the side of it. The air convection/con- 
duction correction terms were measured by radi- 
antly and electrically heating the ECR in air and 
vacuum. These experiments have been described 
by Geist [30], The only correction which was mea- 
surably non-zero was the air convection/conduc- 
tion term. It amounted to 0.1%. 

To estimate the error caused by spatial non-uni- 
formity of response, we scanned the ECR aperture 
with a 1-mm diameter, 514-nm laser beam, and also 
mapped the radiation field produced by the laser 
sphere at this aperture. The former was found to 
vary by 0.6% (see fig. 13), and the latter by 
-fO.01%. The corresponding correction factor was 
evaluated as 0.9998-1-0.0005. 

Since the ECR is presumed to be sensitive in the 
infrared, it was necessary to determine the contri- 
bution which heating of the laser-sphere source 
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Figure 13. Relative horizontal distribution of ECR responsivity 
(514 nm). 

makes to the ECR signal. This contribution was 
estimated by measuring the radiation temperature 
rise of the sphere interior during laser irradiation at 
5 14 nm, using an imaging infrared radiometer. The 
measured temperature rise was 0.09 K. Hence, the 
increment in the radiant power received by the 
ECR was calculated as 0.056 fiW, about 0.01% of 
the 0.4-mW radiant power levels used in our exper- 
iment. 

The results of these various measurements and 
their estimated, three-sigma uncertainties are listed 
in table 4. We applied a 0.1% correction to the 
ECR responsivities measured by electrical substitu- 
tion heating and estimated the overall systematic 
uncertainty of the ECR as 0.2%. 

Table 4. Correction factors for ECR spectral responsivity and 
estimated uncertainties (3cr) 



Electrical lead heating 
Radiant case heating 
Reflection and emission loss 
Convection-conduction loss 
Non-uniformity of response 
Infrared response 

Overall correction 



Factor 


Uncertainty (%) 


1.0000 


0.06 


1.0000 


0.05 


1.0001 


0.10 


1.0010 


0.10 


0.9998 


0.05 


1.0001 


0.01 



1.0010 



0.17 



4.7 Monitor Diode and Neutral Density Filter 

Any measurement of the laser sphere by any of 
the three detector systems was accompanied by a 
simultaneous measurement by the monitor detec- 
tor, followed by dark-current readings obtained by 
blocking the laser input to the sphere. The ratio of 
the dark-current corrected detector signal to the 
dark-current corrected monitor signal was the 



quantity actually used in the subsequent analysis. 
As mentioned earlier, this suppressed errors due to 
fluctuations of laser-sphere output. The monitor 
diode was an EG&G 444B silicon photodiode. Its 
response was found to be stable within 0.01% over 
periods of several hours. 

Because the optimal power for ECR measure- 
ments is about four orders of magnitude higher 
than the power level produced by the blackbody, a 
neutral density filter was inserted in the laser beam 
when silicon-diode and spectroradiometer readings 
were taken. The monitor-diode signals obtained 
with and without the filter provided measurements 
of the transmittance of this filter, so that our data 
were automatically adjusted for its presence when 
ratios were computed. 

4.8 Detector Linearity 

The linearity of response of the spectroradiome- 
ter and silicon detectors was measured with an au- 
tomated beam conjoiner [31], a flux-addition 
device. For the range of signal levels and the oper- 
ating conditions used in these gold-point measure- 
ments, we found no correction necessary for 
departures from strict linearity. The ECR was 
tested by varying the electric heater power and ex- 
hibited no detectable nonlinearity over the small 
dynamic range in which this instrument was used. 
We assessed the combined uncertainty due to de- 
tector nonlinearities to be less than .0.1%. 

4.9 Apertures 

The mutual arrangement of the laser-sphere and 
detector apertures is shown in figure 14. The laser- 
sphere aperture was oriented perpendicular to the 
spectroradiometer optic axis by turning on the 
alignment laser in the spectroradiometer, holding a 
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Figure 14. Laser-sphere and detector apertures. 
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glass microscope slide against the aperture, and ad- 
justing the sphere until the laser beam was re- 
flected by the glass back along itself. In a similar 
manner, the detector apertures were oriented using 
an auxiliary HeNe laser mounted on the source 
table but aimed into the spectroradiometer pre- 
cisely along the spectroradiometer alignment laser 
beam. Because of the precision of these alignments, 
we believe that any error arising from non-paral- 
lelism of the apertures is negligible. These co- 
aligned lasers also provided means of positioning 
the centers of the detector and sphere apertures in 
the same horizontal plane, and of determining the 
spectroradiometer carriage positions at which the 
detectors are opposite the laser-sphere. 

Reflections from the diode surface and from the 
diode and ECR apertures and mounting hardware 
could potentially be returned by the sphere aper- 
ture or hardware back to the detectors, or be scat- 
tered within the sphere and partially returned to 
the detectors. Although numerical estimates sug- 
gested that the errors due to such reflections are 
exceedingly small, two verification experiments 
were carried out. In the first, the detectors were 
covered by black velvet cloth except for the aper- 
tures and a narrow sliver at their edges. No de- 
tectable (<0.01%) signal differences were 
observed. In the second, the monitor signal was 
followed as the detectors were moved into position 
and away. Again, there was no detectable signal 
change. 



5. Theory 

5.1 Simplified Derivation of Measurement 
Equation 

The basic principle of our measurement is best 
explained by temporarily neglecting small effects 
arising from the finite bandwidth of the sources 
and detectors involved, as well as all further effects 
introduced by their different sizes and spatial radi- 
ance distributions. With these simplifications, the 
spectroradiometer signals obtained at any given 
wavelength with the gold-point blackbody and 
laser-sphere sources are given by 



S.f= 



--LrAXGsRs, 



SLS r LS /^ D 



(13) 
(14) 



where Lx° AX and L^, respectively, are the radi- 
ances of the blackbody and the laser sphere, AX is 



the spectral bandwidth of the spectroradiometer, 
and Gs and Rs, respectively, denote the geometric 
extent and the spectral responsivity of the spectro- 
radiometer system. By taking the ratio of eqs (13) 
and (14), we obtain 



Lf AX=(SiVSk^L''^ 



(15) 



thereby eliminating the unknown properties of the 
spectroradiometer and relying only on the linearity 
of its response to ensure that this result is accurate. 
Similar to eq (14), the signal obtained with the sili- 
con-diode or ECR reference detector and the laser 
sphere is 



cLS r LS /^ r> 

Od — JLj Ltd -f^Di 



(16) 



where i?D is the spectral responsivity of the detec- 
tor, and where the geometric extent Go is given by 
the standard expression which describes the radia- 
tion transfer between circular, parallel, and coaxial 
apertures with radii /"ls and /-q, and separation d 
[32], 



GD = (7rV2) 



(d' + rls+rh)-Kd'+rls+,iy 



-^rlA] 



1/2 



(17) 



Combining eqs (15) and (16), we obtain the follow- 
ing approximate result for the spectral radiance of 
the blackbody in terms of the measured spectrora- 
diometer and reference detector signals, the geo- 
metrical parameters of the laser sphere and 
reference detector arrangement, and the known 
spectral responsivity of the reference detector. 



LZ==:[SrSh'/(Sk'A\RuGu)]. 



(18) 



The subscript (0) has been affixed to this expression 
for Lx° in order to indicate that it is only an ap- 
proximation to the full measurement equation. The 
corrections which must be applied to it for the rig- 
orous reduction of measured data analysis will be 
derived next. 

5.2 Derivation of Correction Factors 

5.2.1 Spectroradiometer Signals The rigorous 
measurement equation describing the spectrora- 
diometer signals obtained when the centers of the 
blackbody and laser-sphere source are imaged with 
unit magnification upon the small rectangular en- 
trance slit of the instrument is [33]: 
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•Hix,y)dxdy. (19) 

Here, 

5s(X) is the spectroradiometer signal; 

X is the wavelength setting of the monochromator 
and X' is an arbitrary wavelength; 

X and y are cartesian coordinates in the source 
plane, referred to the source center (0,0); 

L\(\'^,y) is the spectral radiance of the source; 

i?s(X') o-(A.,X') is the spectral response of the spec- 
troradiometer to radiation of wavelength X' when 
the monochromator is set at X; 

H(x,y) is the relative spatial response of the spec- 
troradiometer to radiation from a source point 
(x ,y) when the instrument is focused at the source 
center. 

The integration in eq (19) and later expressions ex- 
tends over all wavelengths and over the entire 
source aperture plane. 

The factorization of the spectral response of the 
instrument into the two factors Rs and cr is 
imiquely defined if we require that RsQ^") is inde- 
pendent of X and that the maximum value of 
cr(X,X') is exactly unity. Then cr(X,X') is what is 
commonly known as the slit function shown in fig- 
ure 7 (or, better, the slit-scattering function to em- 
phasize the fact that the distant wings are due to 
diffraction and scattering), and Rs(X') is the spec- 
tral responsivity of the instrument shown in figure 
9, which is dominated by the spectral transmittance 
of the monochromator and the spectral response of 
the detector. With perfect imagery, the spatial ap- 
paratus fimction H(x,y) would be a constant every- 
where within a 0.6x0.8-mm rectangle and zero 
elsewhere. However, due to diffraction and scat- 
tering from dust and optical imperfections we must 
allow for the departures from this ideal illustrated 
by figure 8. We assumed that, like the slit-scatter- 
ing function, H(x,y) is normalized so that its maxi- 
mum value is one, although this is not necessary 
since any such normalization factors will eventu- 
ally cancel out in our application. 

We replace the blackbody spectral radiance by a 
product of its value at the center and a relative 



spatial distribution function. The former is ex- 
pected to be a Planckian function of wavelength. 
The latter is shown in figure 5. It exhibits a value 
close to unity over a small (1.5-mm diameter) cen- 
tral circular area, and then decreases to near zero at 
a radial distance of about 20 mm. Accordingly, we 
assume 



L,(X>,j;)=Lx^«(X') /BB(x,y), 



(20) 



where /bb(0,0)=1. This factorization of the spec- 
tral radiance into separate spectral and spatial fac- 
tors is an acceptable approximation for the 
purposes of the integrals in eq (19) because the slit- 
scattering and spatial apparatus functions are both 
sharply peaked functions. Thus, the error in this 
approximation at large distances from the source 
center or at wavelengths far from the monochro- 
mator setting will make a negligible contribution to 
the value of the integrals. With this approximation 
to the spectral radiance, we can write the measure- 
ment equation for the spectroradiometer viewing 
the center of the blackbody source as 



Si«(X)= /Z,?«(X') i?s(X') o-(X,X') dX' 
■fffBB(.x,y)H(x,y)(ixdy. 



(21) 



We also assume that the spectral radiance of the 
blackbody is Planckian, so that 



Lf(X')- 






PQ^',T), 



(22) 



where Xl is the wavelength of the laser-sphere 
source, Tis the blackbody temperature, and P(X,T) 
is the Planck function defmed by eq (3). Thus we 
rewrite eq (21) as 

5r(?^)=^^^ /^(X',31 /?s(X') cr(X,X') dX' 



•// ^BB{x,y) H(x,y) dx dy. 



(23) 



In order to estimate the possible error introduced 
by the approximations of eqs (20) and (22), we 
modeled the distributions involved by simple ana- 
lytical functions which permitted the exact evalua- 
tion of the integrals in terms of elementary 
functions. By varying the parameters of these mod- 
els to correspond in varying degrees to the mea- 
sured functions, and observing the effect of the 
product approximation on eq (20), we found that 
the error is on the order of 0.002%. 
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We approximate the laser-sphere radiance in eq 
(19) by 

L^{X'^,y) = S(X'-XOL'-KO,0) /lsC^,;;), (24) 

where the delta function reflects the fact that the 
spectral distribution is that of a sharp laser line at 
wavelength Xl, and L^^(0,0) is the central radiance 
of the source. The remaining factor <f Ls(^>y) is the 
relative spatial radiance distribution of the source 
shown in figure 6, which is approximately unity 
over the circular exit aperture of the sphere and 
zero outside this aperture. By performing the inte- 
gration over wavelength in eq (19), we obtain 



< <^Ls> =1— -JI ^i.s{x,y) dx dy. 



(27) 



■I f^us(x,y) H(x,y) dx dy. 



(25) 



5.2.2 Absolute Detector Signals The rigorous 
measurement equation for the absolute (silicon- 
diode or ECR) detector signals is [33] 

Sh^ = Iin IU(X'^ ,y) Ro(X') dX' 

COS0-COS0' J J J < J < /n/rv 

■ ,2 , , K2~r? K2 dx dy dx dy . (26) 

In this expression, 

^ is the detector signal; 

{x,y) and {x'^y"^ refer to parallel cartesian coordi- 
nate axes in the source and detector aperture 
planes, respectively; 

Q and Q' are the angles between the normals to 
these aperture planes and the ray defined by the 
points (x,j;) and {x' yy, 

d is the distance between the aperture planes; 

Lx(X',;c,j;) is the spectral radiance of the laser- 
sphere source, assumed directionally uniform 
(Lambertian) within the small solid angle sampled 
by the detectors (/'/19 for the silicon detector and 
//25 for the ECR); 

i?D(\') is the spectral responsivity of the detector. 

We assume again that the spectral radiance of 
the laser-sphere source is given by eq (24). We fur- 
ther assume that in the integrals of eq (26) the rela- 
tive spatial distribution function can be replaced by 
its average value. 



and moved outside the integrals, ^ls is the laser- 
sphere aperture area. Thus eq (26) becomes; 

^B'=^d(M i'^'(0,0)< /ls> 



•M: 



d-^i^T-xT^iy-yJ^^y ^'^y'- (28) 



For the geometry we used, this integral can be 
evaluated to give 



5B^=i?D(AL)Z-'-'(0.0) < /ls> C7d, 



(29) 



where Go is the geometric extent previously de- 
fined by eq (17). 

5.2.3 Definition of Correction Factors We can 

now combine eqs (23), (25), and (29) in order to 
derive the rigorous analog of the approximate ex- 
pression (18) for the blackbody spectral radiance. 
This is accomplished by dividing eq (25) into eq 
(23), solving for Lx^(\l). substituting the result for 
L'-^(0,0) obtained from eq (29), and then multiply- 
ing the numerator and denominator by the spectral 
bandwidth A\l of the spectroradiometer at X^^. 
Thus we obtain 

rBB^-X >. iSs Q0> 5d (\l) fa-fw f,ri\ 

^x (Al)-^ls(^) ^y^^ R^(j^^) fj^ ■ ^ ^^^ , (30) 

where the bandwidth LX is defined as the normal- 
ized area of the slit-scattering function. 



o-(A,Xl) 



(31) 



The left-hand factor in eq (30) is identical to the 
approximate radiance of eq (18), and the right-hand 
factor is comprised of three correction terms; 

(1) A spectroradiometer slit-scattering function 
correction, defined as 

f P{X^,T)Rs{X^)!cr{X,X')dX' 

^^ fp(x;r)Rs{X')o-{x,x')dX- ■ ^^"^^ 

This factor arises from the different spectrora- 
diometer response to broadband blackbody radia- 
tion and monochromatic laser light. 

(2) A spatial apparatus function correction, defined 
as 
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/h= 



/ / / Ls(-^ ,y) H{x ,y)dx dy 
/ / / bb(^ ,;') H{x ,y) Ax dy ' 



(33) 



This correction accounts for the different sizes and 
spatial radiance distributions of blackbody and 
laser-sphere sources. 

(3) A sphere non-uniformity correction, defined as 
the inverse of the average relative sphere radiance 
of eq (27). This correction accounts for the fact 
that the spectroradiometer viewed primarily the 
center of the sphere aperture, while the silicon and 
ECR detectors viewed the entire aperture. 

5.2.4 Diffraction Effects The theory pre- 
sented above does not include all diffraction effects 
which constitute a potential source of error in the 
physical model leading to eq (30). For the spectro- 
radiometer measurements, diffraction is implicit in 
the slit-scattering and spatial apparatus functions 
considered. Diffraction losses at the limiting aper- 
ture of the spectroradiometer were eliminated by 
limiting the field of view, not at the laser sphere or 
blackbody, but within the spectroradiometer sys- 
tem, so that the effects would cancel when the ra- 
tio 5sV5P in eq (31) is taken. However, for the 
silicon-detector and ECR measurements, there may 
be a small diffraction loss caused by the apertures 
being mounted a few millimeters in front of the 
detectors themselves. We performed calculations 
based on the treatment of Steel, De, and Bell [34] 
and found losses ranging from 0.003% for the sili- 
con detector at 514 nm to 0.02% for the ECR at 
647 nm. The upper limit of this range (0.02%) was 
included in our error analysis. 



6. Measurements 

6.1 Spectroradiometer Slit-Scattering Function 
Correction, f^ 

This correction factor was calculated from eq 
(32) by combining the measured values of the spec- 
troradiometer spectral responsivity RsQ^) (fig- 9) 
and slit-scattering function cr(A.,A.') (fig. 7) with cal- 
culated values of the Planck function P(k,T) [eq 
(3)] evaluated at the IPTS-68 gold-point tempera- 
ture, 1337.58 K. For these calculations we used the 
0.6-mm slit-scattering functions obtained in our 
earlier study, because they include far-wing data. 
Smce the slit-scattering function appears in inte- 
grals in both the numerator and denominator of f^.. 



it tends to cancel and small differences between the 
functions obtained previously and the current slit- 
scattering functions are of no consequence. The 
values of /o- at the three laser wavelengths fall 
within the range 0.9999+0.0004. 

6.2 Spectroradiometer Spatial Apparatus Function 
Correction, /h 

To compute the value of /h defined by eq (33), 
we used the normalized horizontal and vertical tra- 
verses of the spatial apparatus function of the spec- 
troradiometer shown in figure 8 and the measured 
spatial radiance distributions of the blackbody and 
laser-sphere sources shown in figures 5 and 6. 
These distributions, of course, represent convolu- 
tions of the true distributions with the spectrora- 
diometer point-spread function, whereas the 
functions needed to calculate fa are the true distri- 
butions. In the evaluation of the integral in the de- 
nominator of eq (33) only the central 1 or 2 mm 
contribute significantly, and here the effect of con- 
volution is merely a reduction of the apparent radi- 
ance by a nearly constant factor (e.g., 0.9997). 
Consequently, the normalized observed distribu- 
tion is virtually indistinguishable from the true dis- 
tribution insofar as the evaluation of the 
denominator in eq (34) is concerned, and no decon- 
volution is required. A similar remark applies to 
the laser-sphere integral appearing in the numera- 
tor, although in this case no deconvolution correc- 
tion would have been expected anyway due to the 
uniformity of that distribution. The values obtained 
for/n range from a minimum of 1.0002 at the end of 
a melt at any wavelength to a maximum of 1.0008 
at the end of a freeze at 514 nm. Our final result for 
/h is the average of all calculated values, 
/h= 1.0003+0.0006, where the uncertainty covers 
the dependence upon time and wavelength and in- 
cludes an uncertainty for the contribution of the far 
wings of the spatial apparatus function, extrapo- 
lated beyond ±7 mm. 

6.3 Laser-Sphere Nonuniformity Correction, 

</ls> 

The integration of the measured distribution 
over the laser-sphere source (fig. 6) was performed 
at all three wavelengths used. The average relative 
radiance obtained was < /ls> = 1.0022+0.0009, 
where the uncertainty is a three-sigma estimate of 
repeatability. 
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6.4 Geometric Extent of Absolute Detectors, Go 

For the evaluation of this quantity, eq (17) was 
used in a simplified form in which the areas of the 
laser-sphere and detector apertures were used in- 
stead of their radii. The aperture dimensions were 
measured on a movable stage under a microscope 
calibrated by the Precision Engineering Division of 
NIST. The microscope was focused on the edge of 
the aperture being measured, and then the (x,y) co- 
ordinates of the stage were recorded as the micro- 
scope cross-hairs were positioned at arbitrary 
intervals of about 5° or 10° around the edge of the 
aperture. The data points recorded in this way 
were sufficiently numerous (about 50) that the 
omission of half of them by using only alternate 
points resulted in a totally negligible difference in 
the area determination. The (x,y) coordinate pairs 
were converted to {r,Q) coordinate pairs about an 
origin near the aperture center (chosen so that all r 
values were equal within 0.01 mm), and the areas 
were calculated from 



'=J^" 



A=\ rMe/2=2^'(Q'+'-Q'->)/4' 



(34) 



with Qo=On and 9„+i=9i. These aperture mea- 
surements were repeated two or three times for 
each aperture by each of two operators with excel- 
lent agreement. The values obtained and their re- 
peatabilities (three standard deviations) are given in 
table 5. 

Table 5. Aperture dimensions and distances 



Laser-sphere aperture area 


(2.5523+0.0013) cm^ 


Laser-sphere aperture thickness 


(0.015+0.001) cm 


Si-detector aperture area 


(0.50021 ±0.00035) cm^ 


ECR aperture area 


(0.2835±0.0004) cm 


Aperture distance 


(15.000±0.002) cm 


Geometric Extent, Gsi 


(5.63843 ±0.00052) 10"^ cm^sr 


Geometric Extent, Gecr 


(3. 19662 ±0.00049) 10"^ cm^sr 



The separation between source and detector 
apertures was set each morning by holding a paral- 
lel-faced gauge block of known length (150.0 mm 
for most measurements) between the apertures, and 
adjusting the detector position until both apertures 
were in light contact with the gauge. In calculating 
the geometric extents of the two detector systems, 
we accounted for the fact that the laser-sphere 
aperture had a thickness of 0.15 mm and was bev- 
eled toward the inside of the sphere by using an 
effective aperture separation of 150.15 mm. The 



uncertainty in the length of the gauge block was 
assessed as 0.02 mm, and the possible error due to 
lack of sharpness of the aperture edges was esti- 
mated as 0.01 mm. 

6.5 Spectral Radiance Measurements 

The measurement of a single blackbody spectral 
radiance value constituted a day's work which was 
carried out in accordance with the measurement 
protocol described below. The measurement was 
performed in a time-symmetrical fashion, so that 
possible drifts of the spectroradiometer blackbody 
signals would be eliminated in part when overall 
averages were taken. However, no drifts were ob- 
served outside of statistical fluctuations. All laser- 
sphere readings were recorded relative to 
simultaneously taken monitor-diode readings, so 
that errors due to small fluctuations of laser -sphere 
output were minimized and the transmittance of 
the neutral-density filter was automatically ac- 
counted for when it was inserted in the laser beam 
entering the sphere. 

The individual measurement steps performed were: 

(1) While the blackbody was heating up to the first 
melt of the day, the laser was turned on and the 
separation d between the laser sphere and detector 
apertures was set as described in section 6.4. The 
ECR was positioned in front of the laser sphere, 
and the laser power was adjusted to give readings 
of approximately 0.4 mW (without neutral density 
filter). 

(2) The ECR was calibrated by applying electrical 
substitution heating for a duration of 180 s, and 
then its responsivity R^cr (V/W) was recorded as 
the ratio of the thermocouple EMF and heating 
power applied, divided by the ECR correction 
1.001. 

(3) The ECR was exposed to radiant heating for 
180 s, and the relative ECR laser-sphere signal S^cr 
was recorded (in V/A) as the ratio of the thermo- 
couple EMF and the monitor-diode reading. 

(4) The neutral-density filter was inserted in the 
beam and the spectroradiometer positioned in front 
of the laser sphere. The central 10-nm portion of 
the spectroradiometer slit-scattering function was 
recorded by scanning the monochromator wave- 
length drive (5 14 and 633 nm) or scanning the dye 
laser (647 nm). The monochromator was set to the 
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maximum of this scan, where it remained until the 
measurements were complete for the day. A pre- 
liminary value, AXo. of the spectroradiometer band- 
width was calculated from eq (31) by integrating 
the measured slit-scattering function. This result 
was corrected for the omitted contribution of the 
far wings of the slit-scattering function as described 
in section 4.4. 

(5) Once the blackbody was in a melt (or, later, in a 
freeze), a cycle of three readings was taken repeat- 
edly: the spectroradiometer readings of the black- 
body and the laser sphere, and the silicon and 
monitor diode readings of the sphere. The spectro- 
radiometer blackbody signal 51^ was recorded as 
the measured photomultiplier anode current (A). 
The spectroradiometer and silicon-diode laser- 
sphere signals, 5s^ and Ssi^, were recorded relative 
to the monitor-diode readings (in dimensionless 
units). 

(6) During the hour or so following a melt (or 
freeze) and preceding the subsequent freeze (or 
melt), the laser sphere was repeatedly measured 
with the ECR and the silicon and monitor detec- 
tors, as described in steps (3) and (5). Generally, 
two to four melt-freeze cyles were carried out in a 
day, each time repeating the measurements of step 
(5). 

(7) Following the last freeze of the day, the laser 
sphere was re-measured with the ECR and the sili- 
con and monitor detectors, the spectroradiometer 
slit-scattering function was re-measured as in step 
(4), and the ECR was re-calibrated by electrical 
substitution as in step (2). 

Table 6 gives a representative sample of the av- 
erage data obtained at 647 nm on a typical day. The 
table also gives the values of the various other fac- 
tors needed to calculate the blackbody spectral ra- 
diance from eq (30). In the example shown, the 
silicon-detector and ECR results are 63.319 and 
63.395 Wcm'^sr"', respectively. 

6.6 Blackbody Spectral Radiances 

The measurement sequence described in section 
6.5 was performed five times at each laser wave- 
length. This gave 10 spectral radiance values, each, 
at the argon and krypton laser wavelengths, and 
five values at the helium-neon laser wavelength. 
The 25 spectral radiance values thus obtained are 
listed in table 7 and were treated as single measure- 
ments for data reduction purposes. 



Because our experiment relies on the redundant 
detector scales established by the silicon diode and 
the ECR, we verified the mutual agreement of 
these scales by comparing the differences and the 
standard deviations of the silicon-diode and ECR 
values of table 7. As shown at the bottom of the 
table, the two scales were found to agree within 
experimental uncertainties. 

Table 6. Representative data set obtained at 647.1 nm 



ECR responsivity, Recr 

Relative ECR signal, SkcR 

Spectroradiometer bandwidth 
(including far wings), AX 

Spectroradiometer blackbody 
signal, SP 

Relative spectroradiometer 
laser-sphere signal, Sp 

Relative silicon-diode signal, SsT^ 

Silicon-diode responsivity, Rs\ 

Silicon-diode geometric extent, Gsi 

ECR geometric extent, Gecr 

Spectroradiometer slit-scattering 
correction, fcr 

Spectroradiometer apparatus-function 
correction, fn 

Laser-sphere nonuniformity 
correction, < /ls> 

Blackbody spectral radiance (Si), L™ 

Blackbody spectral radiance (ECR), Lj?^ 



0.267244 V/W 
0.19412 V/W 
2.5273X10-'cm 

3.17858X10-* A 

0.449073 

0.599112 
0.468170 A/W 
5.63843 X 10-^ cm^sr 
3.19662xlO-'cm^sr 
0.9999 

1.0003 

1.0022 

63.3186 Wcm-^sr-' 
63.3953 Wcm-'sr-' 



Table 7. Measured blackbody spectral radiances (Wcm~'sr~') 
and random uncertainties (Itr) 



514.5 nm 632.9 nm 
Si ECR Si Si 



647.1 



ECR 



Mean 



2.7468 2.7523 48.883 63.482 63.484 

2.7576 2.7570 48.621 63.319 63.395 

2.7656 2.7612 48.761 63.192 63.178 

2.7650 2.7605 48.659 63.396 63.379 

2.7653 2.7616 48.696 63.477 63.231 

2.7601 2.7586 48.724 63.373 63.333 



Stddevofmean 0.0036 0.0018 0.046 0.054 0.056 
(0.13%) (0.07%) (0.09%) (0.09%) (0.09%) 



Difference 



0.0015 
(0.05%) 



0.040 
(0.06%) 
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The standard deviations of the mean listed in 
table 7 represent the composite random uncertain- 
ties of all of the measurements described in section 
6.5. They include the random uncertainties of all 
detector signals, the calibration of the ECR, the 
wavelength setting of the spectforadiometer and 
measurement of its bandwidth, the setting of the 
aperture distance, the transmittance of the neutral 
density filter, as well as fluctuations of the black- 
body and laser-sphere sources. Since the individual 
standard deviations are all on the order of 0.1%, 
our three-sigma estimate of the composite random 
uncertainty of these measurements is 0.3%. 

6.1 Gold-Point Temperature 

To find the gold-point temperature Tau, we ex- 
pressed the theoretical blackbody radiance as the 
product of Planck's eq (3) and the estimated black- 
body emissivity, e= 0.9999, of section 4.1. We com- 
pared the average measured spectral radiances Lf^ 
given in table 7 with the values e P(X,T/^a) obtained 
by substituting trial values of T^u, the laser wave- 
lengths listed in table 3, the refractive indexes of air 
computed from Edlen's formula [35], and evaluated 
the radiation constants, C\ = 2hc^ and C2=hc/k, in 
terms of the 1986 values of A, k, and c given by 
Cohen and Taylor [36], By varying Tau until the 
quadrature sum of the weighted differences 
[Z-K^-e P(X,rAu)]/o-\ (o-K being the standard devia- 
tions of the mean given in table 7) was minimized, 
we found Tau^ 1337.334 K. The relative standard 
deviation of the mean of this least-squares fit vvas 
±0.056% (relative to radiance). The theoretical 
blackbody spectral radiances eP(X,T/^a) calculated 
for this temperature and their deviations from the 
average measured spectral radiances Lf^ are listed 
in table 8. These residuals exhibit a small depen- 
dence on wavelength which is unexplained, but 
falls within the combined uncertainty of our mea- 
surements. 

Table 8. Residuals of least-squares fit for rA„= 1337.334 K 

Wavelength Calculated radiance Residuals 

(nm) (Wcm-'sr-') Si ECR 



514.5 
632.9 
647.1 



2.7561 
48.736 
63.433 



0.0040 
-0.12 
-0.060 



0.0025 



-0.100 



6.8 Uncertainty 

To estimate the uncertainty of this result, we 
added in quadrature the composite random uncer- 
tainty of the blackbody spectral radiance measure- 
ments (sec. 6.6), the uncertainty of the least-squares 
fit (sec. 6.7), and the various systematic uncertain- 
ties in excess of 0.01% that have been mentioned in 
previous sections of this paper. These uncertainties 
have been collected in table 9, where all of them 
are given as estimated upper limits of systematic 
errors or three standard deviations. Their summa- 
tion in quadrature gives a total uncertainty (relative 
to radiance) of 0.45%. Hence, the uncertainty of 
the measured gold-point temperature was esti- 
mated from eq (10) as ±0.34 K (at wavelength 600 
nm). 

Table 9. Estimated uncertainties (3cr) of gold-point measure- 
ment (see sections cited in parentheses for details) 



Quality 



Uncertainty 



Composite random uncertainty (6.6) 0.30% 

Least-squares fit (6.7) 0.17% 

Blackbody emissivity (4.2) 0.01% 

Polarization (4.3) 0.02% 

Far-wing bandwidth correction (4.4) 0.04% 

Detector responsivity [silicon (4.5) or ECR (4.6)] 0.20% 

Detector linearity (4.8) 0.10% 

Diffraction (5.2.4) 0.02% 

Slit-scattering function correction /^ (6.1) 0.06% 

Apparatus function correction /h (6.2) 0.03% 

Sphere non-uniformity correction </ls> (6.3) 0.09% 

Geometric extent, Go (6.4) 0. 15% 

Total radiance uncertainty (added in quadrature) 0.45% 

Gold-point temperature uncertainty 0.34 K 
[from eq (10), at 600 nm] 



7. Result 



Based on the above, the value 
rAu=(1337.33±0.34) K 



(35) 



is quoted as the final result of our measurement of 
the gold-point temperature.^ This result is in good 
quantitative agreement with the measurements 



^This value differs from a preliminary value which was re- 
ported at an IMEKO Photon-Detector Committee meeting 
(Braunschweig 1987). The preliminary value was based on an 
incomplete analysis of the raw data reported in the present pa- 
per, and was in error because the refractive index of air (about 
1.00027) had not been included in Planck's radiation formula. 
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described in section 3, above. Our experiment rep- 
resents the only direct spectroradiometric determi- 
nation of the gold-point temperature which has 
been reported so far and, therefore, provides an 
independent confirmation of these earlier measure- 
ments. The value quoted in eq (34) was submitted 
to the Consultative Commmittee on Thermometry 
as an official NIST contribution to the ITS-90 and, 
as noted in section 1, was incorporated into the 
new scale. 
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